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Abstract 

We report observations of radiative B meson decays into the K + ir~~f and K + ir~~ir + j final states. 
In the B° — > K + ir~j channel, we present evidence for decays via an intermediate tensor meson 
state with a branching fraction of B(B° -> K|(1430)°7) = (1.3 ± 0.5(stat.) ± 0.1(syst.)) x 1CT 5 . 
We measure the branching fraction B(B + — > i"C + 7r~7r + 7 ) = (2.4 ± 0.5(stat.) l^(syst.)) x HT 5 , 
in which the i? + — > i^T* 7r + 7 and -B + — > K + p°~/ channels dominate. The analysis is based on a 
dataset of 29.4 fb _1 recorded by the Belle experiment at the KEKB collider. 

PACS numbers: 13.40.Hq, 14.40.Nd 
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Since the first measurement of the inclusive branching fraction for B — > Xg^y by the CLEO 
collaboration in 1995 0, the flavor changing neutral current process b — > s'j has been used 
as a sensitive probe to search for physics beyond the Standard Model (SM). In experiments 
at the T(AS), a pseudo-reconstruction technique, in which the X s state is reconstructed from 
one kaon and multiple pions, has been the most powerful tool to identify b — > 57 events. 
In order to measure more precisely the inclusive rate, a detailed knowledge of the exclusive 
final states is required. In addition to the already established B — > K* r y decay 0, there 
are several known resonances that can contribute to the final state. CLEO has reported 
evidence for B — ► K\ (1430)7 [0]- Some theoretical predictions for the branching fractions 
of the exclusive decays can be found in Ref. [|J. Exclusive decays, such as B — > K\ (1400)7, 
can also be used to measure the photon helicity, which may differ from the SM prediction 
in some new physics models ||. 

In this Letter, we report on a search for resonant structures Kx above the K* mass in 
radiative B meson decays. The analysis is based on a data sample of 29.4 fb _1 (31.9 million 
BB events) recorded by the Belle detector || at KEKB . KEKB is an asymmetric energy 
e + e _ collider (3.5 GeV on 8 GeV) operated at the T(45) resonance. The Belle detector 
has a three-layer silicon vertex detector (SVD), 50-layer central drift chamber (CDC), an 
array of aerogel Cherenkov counters (ACC), time-of-flight scintillation counters (TOF), an 
electromagnetic calorimeter of CsI(Tl) crystals (ECL). 

We select events that contain a high energy photon (7) with an energy between 1.8 and 
3.4 GeV in the T(4S) center-of-mass (CM) frame and within the acceptance of the barrel 
ECL (33° < # 7 < 128°). In order to reduce the background from Tr°,r] — > 77 decays, we 
combine the photon candidate with all other photon clusters in the event and reject the 
candidate if the invariant mass of any pair is within 18 MeV/c 2 (32 MeV/c 2 ) of the nominal 
7T° (rj) mass (this condition is referred to as the 7T°/r] veto). 

We search for Kx resonances decaying into two-body (K + n~) and three-body (K + 7t~tt + ) 
final states pi in the invariant mass (Mk x ) range up to 2.4 GeV/c 2 . For the K + tt~ final 
state, the range Mk x < 1.2 GeV/c 2 is excluded to remove K* contributions. Charged tracks 
are required to have CM momenta greater than 200 MeV/c, and to have impact parameters 
within ±5 cm of the interaction point along the positron beam axis and within 0.5 cm in 
the transverse plane. To identify kaon and pion candidates, we use a likelihood ratio that 
is calculated by combining information from the ACC, TOF, and dE/dx (CDC) systems. 
We apply a tight selection with an efficiency (pion misidentification rate) of 83% (8%) for 
charged kaon candidates and a loose selection with an efficiency (kaon misidentification rate) 
of 97% (28%) for charged pion candidates. 

We reconstruct B meson candidates from a photon and a Kx system by form- 
ing two independent kinematic variables: the beam constrained mass Mt, c = 

\l(Kc,Jc 2 ) 2 ~ (\p* Kx +P;\/cf and AE = E* Kx + E* - E^ m , where E^ m is the beam 
energy, and p*, E* Pk x , E* Kx are the momenta and energies of the photon and the Kx 
system, respectively, calculated in the CM frame. In order to improve the Mb c resolution, 
the photon momentum is rescaled so that \p*\ = (E^ eiim — E* Kx )/c is satisfied. 

The largest source of background originates from continuum qq (q = u, d, s, c) production. 
To suppress this background, we use a Fisher discriminant formed from six modified Fox- 
Wolfram moments |I(J and the cosine of the B meson flight direction (cos#^). The moments 
are calculated in the rest frame of the B candidate to avoid a correlation with M^ c |]TlJ . Signal 
and background events are classified according to a likelihood ratio LR = £ sig /(£ si g + £b g ), 
where the likelihood £ S i g (£b g ) is the product of the probability density functions (PDF) 
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FIG. 1: (a) Mk-k (b) M^ c and (c) Icos^hcil distributions for B° — > K + tt~^ candidates. The 
unbinned ML fit results are shown in (a) and (c). The qq backgrounds are subtracted in (c). 
M bc > 5.27 GeV/c 2 is applied in (a) and (c), and 1.25 GeV/c 2 < Mk-k < 1-6 GeV/c 2 is applied in 
(b) and (c). In (a), AE sideband data is scaled to the unbinned ML fit result and overlaid. 



of the Fisher discriminant and cos 9% for signal (background). The PDFs for the Fisher 
discriminant are determined from Monte Carlo (MC) simulations. For cos 9%, we assume a 
1 — cos 2 9* B behaviour for signal events and a flat distribution for continuum background. 
The selection criteria on the likelihood ratio are chosen so that S/ \J S + N is maximized, 
where S and N are (MC) signal and background yields, respectively. The optimized criteria 
retain 68% of the B° — > K + n~'j signal and 42% of the B + — > K + n~n + ^ signal. 

The B decay signal is separated from background, first by applying a requirement on AE 
and then by fitting the M^ c spectrum. If we find multiple candidates with \AE\ < 0.5 GeV 
and M bc > 5.2 GeV/c 2 in the same event, we take the candidate which gives the highest 
confidence level when we fit the Kx decay vertex (best candidate selection). We then select 
candidates with —100 MeV < AE < 75 MeV, which removes 19% and 3% of signal on the 
lower and higher sides, respectively. We define a AE sideband to be 100 MeV < AE < 
500 MeV at M^ c > 5.2 GeV/c 2 , in which we expect negligible signal contribution. 

In the B° — > K + Tt~ r y analysis, we obtain the Mxn distribution shown in Fig. |l|(a). We 



observe an excess around Mk-k — 1-4 GeV /c | i~2|| . The M^ c distribution with 1.25 GeV/c < 



Mk-k < 1-6 GeV/c 2 is shown in Fig. |l|(b). We fit the M^ c distribution to extract the signal 
yield. The distribution for the qq background is modeled by an ARGUS function in 
which the shape is determined from the AE data sideband. The distribution for the signal 
component is modeled by a Gaussian determined from signal MC calibrated by B~ — ► D tt~ 
data. The signal yield is found to be 27 jlf (stat.) ^(syst.) with a statistical significance of 



5. Oct. Here, the significance is defined as y — 21n(£(0)/£ max ), where C max is the maximum 
of the likelihood and £(0) is the likelihood for zero signal yield. 

The observed signal may be explained as a mixture of three components: B° — > 
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fq(1430)°7, B° -> ir*(1410)°7 and non-resonant (N.R.) B° -> ir+vr"7. In order to separate 
these components, we apply an unbinned maximum likelihood (ML) fit to M^ c , the cosine of 
the decay helicity angle (cos^hei) and Mk-k- The expected cos^hei distributions are sin 2 2#hci? 
sin 2 ^h e i and uniform for these three components, respectively. The PDFs for cos #hei and 
Mk-k are determined from the AE sideband data for qq background, from the corresponding 
MC samples for resonant components, and from an inclusive b — > 37 MC sample [IT]] for 
the non-resonant component. The cos#hei PDFs for signals are distorted up to 20% due to 
a non-uniform efficiency The validity of the method is tested with B~ — > D°tt~ data and 
MC. 

The fit results for Mk-k and cos^hei are overlaid in Figs. 0(a) and |l|(c), and summa- 
rized in Table |. We find evidence for radiative decays via an intermediate tensor state, 
B° — ► K% (1430)°7. The f^*(1410)°7 and non-resonant components are not significant, so 
we set upper limits. The 90% confidence level upper limit iV is calculated from the relation 
J"£(n)dn = 0.9 / °° C(n)dn, where £( n) is the maximum likelihood with the signal yield 
fixed at n. 

We estimate the systematic error due to the fitting procedure as follows. For the signal 
shapes in the Mb c and Mk-k distributions, we vary the mean and width parameters in the 
fit within their experimental errors. We also test the validity of the background PDFs by 
replacing them with those obtained from a qq MC sample. We assign the largest deviation 
in these tests as the systematic error of the signal yield. 

The event selection efficiency for B° — > fQ(1430)°7 is (5-0 ± 0.3)% including the sub- 
decay branching fractions. The error includes contributions from photon detection (2.8%), 
tracking (2.3% per track), kaon identification (0.6%), pion identification (0.5%), event se- 
lection including likelihood ratio, 7r / 77 veto and best candidate selection (2.0%) and un- 
certainty of the sub-decay branching fractions (2.4%). Assuming an equal production 
rate for B°B° and B + B~, this leads to a branching fraction of B° — ► K^l^O) ^ of 
(1.3 ± 0.5(stat.) ± 0.1(syst.)) x 10" 5 . 

The result agrees with the predictions based on a relativistic form factor calculation [|J . 
Our result is also consistent with the CLEO measurement |§ when we neglect the non- 
resonant component and assume as they did that the K* (1410)7 component is negligible. 

In the B + —>■ K + -K^-K +r y analysis, we find additional background sources from a MC study. 
Cross feed from B — ► K*j to B + — > K + 7i^7i +/ ~f becomes negligible after removing positively 
identified B — > Kir^j events. The size of the cross feed from other b — > 57 decays, especially 
from those with a 7r° in the final state, is estimated by using the inclusive b — > 57 MC 
sample. The contribution from the b — » c background is estimated by using a corresponding 
MC sample. 

To extract the signal yield, we fit the Mt, c distribution shown in Fig. @(a). In addition to 
a Gaussian and an ARGUS function to describe the signal and qq background components 
obtained using the same method as in the B — > Ktc-/ analysis, smoothed MC histograms for 
the b — > 57 cross feed and other B meson decays are used to model the M bc shape, where 
the normalizations are fixed assuming the luminosity and the measured b — > S7 branching 
fraction [fTT] , TB\. We find the signal yield of 57^ 2 (stat.) ^(syst.) with a 5.9a statistical 



significance. 

The Mk x distribution is shown in Fig. 0(b), where the distribution for qq is obtained 
from the AE sideband and is normalized using the fit result. We observe no signal excess 
above 1.8 GeV/c 2 . The B + — > K + 7T~7T +/ -y signal may be explained as a sum of decays 
through kaonic resonances such as B + — ► Ki(lAQ§) +r y and B + — ► f^*(1680) + 7. The current 
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FIG. 2: (a) M bc , (b) M Kx , (c) M KlT and (d) M W7T distributions. The fit result of the M bc 
distribution is shown in (a), while the result of the unbinned ML fit is shown in (c) and (d). 
M bc > 5.27 GeV/c 2 is applied in (b), (c) and (d). 



statistics and the existence of a large number of resonances prevent us from decomposing 
the resonant substructure. However, it is still possible to measure the K*Tc-f and Rpj 
components separately, as most of the resonances have sizable decay rates through the K*7r 
and Kp channels. 

To find the composition of the signal, we perform an unbinned ML fit to M^ c , Mx n 
and M nn with three signal components (K*TC-f, Kp'j and non-resonant KTiTC-f) and a qq 
background component. In addition, the components from b — > 57 cross feed and from 
other B meson decays are included in the fit with fixed normalizations. The Mk-k and M wn 
shapes for the qq background are determined from the AE sideband data, and those for the 
other components are determined from the corresponding MC samples. 

In order to model the signal PDF for the K*-K / ~f component, we use a mixture of B + — > 
iTi(1400) + 7 -> K*°n + -f and B + -> A^(1680)+7 -> K*°tt + j MC. The K x (1400)7 fraction of 
the mixture is determined to be 0.74 ± 0.14 by examining a background-subtracted M/f m 
distribution for candidates with \M^ — M^* \ < 75 MeV/c 2 {K* mass cut). Likewise for the 
Kp-f PDF, a mixture of B + -> /^(mo)^ -> K+p°7 and B + -> A*(1680) + 7 -> K + p° 7 
MC is used, where the (1270)7 fraction is determined to be 0.68 ± 0.17 according to a 
background-subtracted Mk-k-k distribution for candidates with \M V7T — M p \ < 250 MeV/c 2 
and \M K7T - M K , \ > 125 MeV/c 2 (p mass cut). 

Figures |](c) and ||(d) show the distributions and fit results for Mk-k and M nn . The 
selection efficiency is estimated from a MC sample with the mixture of resonances used for 



the PDF determination. We also consider other well-established resonances |14| which give 
slightly different efficiencies, and assign the difference in the result as a systematic error. 
The signal yields, efficiencies and the branching fractions are listed in Table |. The total 
B + — > i^ + 7r~7T + 7 branching fraction is dominated by B + — > K*°n + ^ and B + -> K+p ^; the 
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TABLE I: Measured signal yields, statistical significances, reconstruction efficiencies, branching 
fractions (B) and 90% confidence level upper limits (UL) including systematic errors. The first 
and second errors are statistical and systematic, respectively. Efficiencies include the sub-decay 
branching fractions Jl5| ]. Efficiencies for K + tt~^/ and K + -K^ir + ^ are based on a mixture of the 
measured sub-components. 



Mode 


Signal Yield 


UL (yield) Sif 


niificance Efficiency(%) 


B (xlO & ) 


UL (xlO 


K + TT~J f 


97 +8 +1 
z ' —7-3 




5.0 § 


18 ±2 


U -^ D -0.12 -0.07 




Kl (1430)°7 


91 +8 +0 

zi _ 7 _ 1 




3.2 


5.0 ±0.3 


1.3 ±0.5 ±0.1 




K*(1410)°7 


7 7+7.I +0.5 
'■' -5.7-1.3 


19 




0.58 ±0.12 




13 


K+TT-j (N.R.) f 


O.Ol^rbO.O 


15 




19 ± 1 




0.26 


K + 7r~7r + 7 | 


cr 7 +12 +b 
°' -11 -2 




5.9 § 


7.5 ±0.7 


2.4 ±0.5^ 




K+p -/ ± 


33 ±^ =b 2 
24 ± 12+^ 


43 


3.7 
2.2 


5.0 ±0.5 
7.4 ±0.7 


2.0t^±0.2 
1.0±0.5i° ; 2 3 


2.0 


K+TT-1T+"/ (N.R.) i 




20 




7.6 ±0.7 




0.92 


Ki(1270)+7 


4.0 ±2.4 ±0.6 


10 




0.40 ± 0.08 




9.9 


ifi(1400)+7 


26±6±o 


36 




2.6 ±0.3 




5.0 



t 1.25 GeV/c 2 < M Kn < 1.6 GeV/c 2 i M Kt7T < 2.4 GeV/c 2 § M bc fit result 



statistical significance for the sum of the two is calculated to be 6.2a and the non-resonant 
component is consistent with zero. We find evidence for the decay B + — > K*°-K +r y with a 
3.7<r significance, while the B + — > K + p 0, y channel alone yields only 2.2a. Systematic errors 
are evaluated using the same procedures as in the B — > Ktcj analysis. 

We also search for resonant decays by applying further kinematical requirements. We 
search for B + — > f^i(1270) + 7 in the K + p 0/ y final state by applying the p mass cut and 
\Mk x — M Kl (i27o) I < 100 MeV/c 2 . We find 6 candidates with a background expectation of 
2.0 ± 0.6 events. To find B + -> iC 1 (1400) + 7 in the K*°7T + j final state, we apply the K* 
mass cut and \M Kx — M^ 1 ( 140 o)| < 200 MeV/c 2 . We obtain a sizable signal; however we only 
provide upper limits due to a lack of ability to distinguish these resonances. The results are 
also listed in Table Q. 

In conclusion, we have studied radiative B decays with the K + n~'y and i^ + 7r _ 7r + 7 final 
states. For K + 7r~j, we consider B° — > K^liSO) ^, B° — > if*(1410)°7 and non-resonant 
components, and find that only the first one is significant. For B + — > K + 7r~7T + ^f, we observe 
the decay mode and measure the branching fraction. The branching fractions for B — > K*n^i 
and Kp'j are consistent with the sum of predicted rates of resonant decays f|] . As listed in 
Table |T|, we find (35 ±8)% of the total B — > X s ^ decay is accounted for by the B — > K*j, 
B -> K* (1430)7, and B -> Kinr'y final states. 
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TABLE II: Exclusive and inclusive branching fractions for the b — > sj process. Equal branching 
fractions are assumed for neutral and charged B decays. Using isospin, the branching fraction of 
B+ -> K* + tt°j (K°p + >y) is assumed to be half (twice) of that for B + if*°7r+7 (K+p ^). 

Mode B (xl(T 5 ) Ref. 



B^K*j 4.2 ±0.4 |,|7| 

B if* (1430)7 (excluding K*ttj, Kpy) 0.9 ± 0.3 

B ^ K*itj 3.1 ±1.0 

B^Kp-i 3.0 ±1.6 

Sum of exclusive modes 11.2 ± 2.1 



B -» X s7 (inclusive) 32.2 ±4.0 0, @ 



Education, Science and Sport of the Republic of Slovenia; the National Science Council and 
the Ministry of Education of Taiwan; and the U.S. Department of Energy. 
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